The Annals of Probability 
2008, Vol. 36, No. 5, 1777-1789 
DOI: 10.1214/07-AOP376 

© Institute of Mathematical Statistics. 2008 

THE LAW OF THE SUPREMUM OF A STABLE LEVY PROCESS 
WITH NO NEGATIVE JUMPS 

By Violetta Bernyk, Robert C. Dalang 1 and Goran Peskir 2 

Ecole Polytechnique Federals, Lausanne, Ecole Polytechnique Federale, 
Lausanne and University of Manchester 

Let X — (X t )t>o be a stable Levy process of index a £ (1, 2) with 
no negative jumps and let St = sup 0<s<t X a denote its running supre- 
mum for t > 0. We show that the density function ft of St can be 
characterized as the unique solution to a weakly singular Volterra 
integral equation of the first kind or, equivalently, as the unique solu- 
tion to a first-order Riemann-Liouville fractional differential equation 
satisfying a boundary condition at zero. This yields an explicit series 
representation for ft. Recalling the familiar relation between St and 
the first entry time t x of X into [a:,oo), this further translates into 
an explicit series representation for the density function of r x . 

1. Introduction. In our study [3] of optimal prediction for a stable Levy 
process X = (X t )t>o, we encountered the question of computing the distri- 
bution function of St = sup 0<s<t X s for t > 0. In the existing literature, such 
expressions seem to be available only when X has no positive jumps and the 
purpose of the present paper is to seek similar expressions when X has no 
negative jumps. We note that the latter problem dates back to [5], page 282. 

Our main result (Theorem 1) characterizes the density function / of S\ 
as the unique solution to a weakly singular Volterra integral equation of the 
first kind or, equivalently, as the unique solution to a first order Riemann- 
Liouville fractional differential equation satisfying a boundary condition at 
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zero. This characterization yields an explicit series representation for / 
(which, in the case of a Brownian motion, coincides with the well-known 
expression arising from the reflection principle). 

Using the scaling property of X, the result extends to St for t 7^ 1. Re- 
calling the familiar relation between St and the first entry time t x of X 
into [x,oo), this further translates into an explicit series representation for 
the density function of t x for x > 0. Moreover, using the Laplace inversion 
formula, we derive an integral representation for / (Corollary 2). Finally, 
we note (Corollary 3) that the proof yields exact constants in the known 
asymptotic expressions for / at zero and infinity. The knowledge of these 
constants plays a key role in our treatment of the optimal prediction problem 
[3]. 

2. The result and proof. 1. Let X = (Xt)t>o be a stable Levy process of 
index a S (1,2) whose characteristic function is given by 

(2.1) _ exp( i f (e - - 1 - a, )fF |^) - .«-«>- 

for A E R and t > 0. It follows that the Laplace transform of X is given by 

(2.2) Ee~ XXt = e txa 

for A > and t > (the left-hand side being +00 for A < 0). From (2.1) and 
(2.2), we see that the characteristic exponent of X equals ^(A) = (— i\) a , the 
Laplace exponent of X equals ip(\) = X a for A > and ip{p) := ip~ 1 (p) = p 1 ^ a 
for p > 0. 

2. The following properties of X are readily deduced from (2.1) and (2.2) 
using standard means (see, e.g., [4] and [13]): the law of (Xct)t>o is the same 
as the law of (c 1 / X t )t>o for each c > given and fixed (scaling property); X 
is a martingale with EXt = for all i > 0; X jumps upward (only) and creeps 
downward [in the sense that P(X Px = x) = 1 for x < 0, where p x = inf{t > 
0:Xt <x} is the first entry time of X into (—00, x)); X has sample paths 
of unbounded variation; X oscillates from —00 to +00 (in the sense that 
liminfj^oo X t = —00 and limsup^^ X t = +00, both a.s.); the starting point 
of X is regular [for both (— 00, 0) and (0, +00)]. Note that the constant 
c = l/r(— a) in the Levy measure v(dx) = (c/x 1+ °) dx of X is chosen/fixed 
for convenience so that X converges in law to y/2B as a f 2, where B is a 
standard Brownian motion, and all the facts below can be extended to the 
general constant c> depending on a if needed (see Remark 2 below). 

3. Let St = sup 0<s<t X s denote the running supremum of X for t > and 
let t x = inf{t >0:Xt>x} be the first entry time of X into [x, 00) for x > 0. 
Since X s > X s _ for all s £ [0,t] and X is right-continuous, one sees that 
P(St > x) = P(t x < t), so the law of St follows from the law of t x (and vice 
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versa). If X is a Levy process with no positive jumps, then it is known that 
the two measures 

(2.3) tP(r x edt)dx = xP(X t edx)dt 

coincide on the Borel cr-algebra of M + x M + (see, e.g., [4], page 190, or [7] and 
the references therein) . This implies that the law of Xt yields the law of t x . It 
follows, in particular, that the known series representations for the density 
function of X t (see, e.g., [17], pages 87-89) lead to series representations 
for the density function of St- If X has no negative jumps, however, then 
the identity (2.3) breaks down and no series representation for the density 
function of St seems to be available in the literature. We mention, however, 
that there is a literature on the distribution of S a when a is an independent 
and exponentially distributed random variable, the process X has arbitrary 
negative jumps, and its positive jumps form a compound Poisson process 
with the jump-size distribution of the so-called "phase type" (see, e.g., [14] 
and [2]). 

4. Our main result can be stated as follows. Note that St == t l / a S\ by 
the scaling property of X so that there is no restriction in assuming that 
t = 1 in the sequel. Recall, also, that ID)" -1 denotes the Riemann-Liouville 
fractional derivative of order a — 1 given by 

(2-4) IT- 1 / (x) = - 1 ~ T 1 f{V } i d y 

for x > and any (admissible) function /:K+ — > M (for more details, see, 
e.g., [16], pages 449-452, and [15], Chapter 2). 



Theorem 1. Let X = (Xt)t>o be a stable Levy process of index a £ (1, 2) 
satisfying (2.1) and (2.2), and let S\ = sup 0<t<1 Xj denote its supremum 
over the time interval [0, 1] . Then the density function f of S\ can be char- 
acterized as the unique solution to the weakly singular Volterra integral equa- 
tion of the first kind 

( 2 -5) f (v + TTTTT" w ^T^r) M ^ 



T(2-a) (x-y) 0-1 / y T(l/a) 
or, equivalently, as the unique solution to the fractional differential equation 

(2.6) xf{x) + aB a ~ 1 f(x) = 
satisfying the boundary condition 

(2.7) limx 2 ~ a f(x) 



xio v ' r(a-l)r(l/a) 
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where D Q_1 denotes the Riemann-Liouville fractional derivative given by 
(2.4) above. This yields the series representation 



00 -I 

(2 - 8) m = £ r( m -i)r(-» t H 1W --'' 

for x > 0. 

Proof. To connect the present result with the existing theory, we will 
begin by recalling a number of known facts about Levy processes with no 
positive jumps (for further details, see, e.g., [4], Chapter VII, and [13], Chap- 
ter 8). 

Let X = (X t )t>o be a Levy process with no positive jumps starting at 
zero, let ^ denote its characteristic exponent, let ip denote the Ijaplace 
exponent of —X and let <J5 := denote the (right) inverse of ip. Thus, 

the characteristic function of X is given by 

(2.9) Ee ai?t = e* J ( A) 

for A £ K. and the Laplace transform of —X is given by 

(2.10) Ee AX t=e t^(A) 
for A > and t > 0. Let 

(2.11) St = sup X s and It = inf X s 

0<s<t 0<s<t 

for t > and set 

(2.12) t x = inf{t > : X t > x} 
for x > 0, on assuming that the infimum is finite a.s. 

From the fact that (e v ^ Xt ~ pt )t>o is a martingale (and that X creeps 
upward), one finds, using the optional sampling theorem, that the Laplace 
transform of t x equals 

(2.13) Ee -PT* = e -xv(p) 

for p > and x > 0. Moreover, if a p is an exponentially distributed random 
variable with parameter p > 0, meaning that P(<r p G dt) =pe~ pt dt for t > 0, 
which, moreover, is independent of X, then (2.13) implies that 

(2.14) P(S ap >x) = P(t x < a p ) = Ee-™ = e~ x ^ 

for p > and x > 0. This shows that S a is exponentially distributed with 
parameter <p(p). Hence, one finds that 

(2.15) Ee A ^- ^ {P) 



(p(p) - A 
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for p > and A £ C with K(A) < p(p). 

Invoking the Wiener-Hopf factorization (see, e.g., [4], page 165, or [13], 
Theorem 6.16) 



(2.16) Ee CTp = Ee p Ee 

V ; p-*(A) 

it follows, using (2.15), that 



(2.17) Ee 



\In-~, _ P(<P(P) ~ A ) 



for A > and p > 0, on recalling that \E'( — iA) = tp(X) for A > 0. The iden- 
tity (2.17) is well known (see, e.g., [4], page 192, or [13], page 213). 

Clearly, X has no negative jumps if and only if X := —X has no positive 
jumps, so, by focusing on the left-hand side of (2.17), one finds 

— /*oo — /*oo 

Ee A/ CTp = p / E ( e A/ t y- P t dt = p E(e- XSt )e- pt dt 
Jo Jo 



(2.18) =p 



1-A / e~ Xx P(S t >x)dx 



~ pt dt 



/■OO rOO 

1-pA/ e~ p '^ / e _Aa: P(5t >x)dx 
Jo Jo 



for A > and p > 0. Combining (2.17) and (2.18) and noticing/recalling 
that ^(A) = = X a and ip(p) = (f(p) =p 1 / Q , one finds that the (joint) 
time-space Laplace transform of (t,x) i— ► P(5< > a;) equals 

/■oo /-cxd 1 /I \a— 1\ 

(2.19) / e~ Xx dx e- pt P(S t >x)dt = — — 7 

V ; Jo Jo p-\ a \p l / a p J 

for A > and p > 0. 

Note that this formula can also be obtained by taking the Laplace trans- 
form with respect to the space variable x on both sides of the expression 



/*oo /*oo 1 

/ e- pt P(S t >x)dt = / e- pt P{T x <t)dt = -Ee 
s Jo Jo p 
(2.20) 

_ ™ p n ~ l x an ^p 



p 



-i F(l + an) , T(an) 

=0 v ' n=l v ' 

where the final identity follows from (8.6) in [13], page 214, combined with 
(ii) and (iii) in [13], page 233. This remark is relevant since the customary 
approach leading to the closed-form expression (2.20) via the so-called scale 
function (cf. [13], pages 214-215) corresponds to Laplace inversion (at least 
formally) with respect to the space parameter. The derivation given below 
takes a different route by firstly performing Laplace inversion with respect 
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to the time parameter and then dealing with the resulting expression using 
techniques of linear integral equations (fractional calculus). 

After these introductory remarks, we are now ready to move to the first 
step of the proof, taking (2.19) as the initial point. Below, we will let L" 1 
denote the inverse Laplace transform with respect to the time parameter 
p and L^" 1 denote the inverse Laplace transform with respect to the space 
parameter A. 

1. Considering p > X a with A > fixed, by (3) in [9], page 238, we find 



(2.211 



l; 



(p - X a )p 1 / c 



(t) 



1 



r(l/a) A 



- 7 (l/a,A a i) 



for t > 0, where (a,x) *— > j(a,x) denotes the incomplete gamma function 



(2.22) 



dy 



for a > and x > 0. Likewise, by (5) in [9], page 229, we find 



A 



a-l 



(2 - 23) ^ L(p-A«)p 

for t > 0. Combining (2.21) and (2.23), we get 



(2.24) L- 1 



1 



1 



A a_1 A] n _l e xat T(l/a,X a t) 
1 A A r(l/a) 



_p - A Q Vp 1/a P 
for t > 0, where (a,x) i— > T(a,x) denotes the incomplete gamma function 



(2.25) 



r(o,x) 



y 



a-l e -y 



dy = T(a) -7(a,x) 



for a > and x > 0. Since the right-hand side of (2.24) defines a bounded 
function of t > and the argument of L" 1 on the left-hand side is a Laplace 
transform defined for all p > [recall (2.19) above], we see that the identity 
(2.24) holds globally for t > and A > 0. 
2. Note that 



1 



(2.26) 



e A t T{l/a,X a t) _ 
X r(l/a) " r(l/a) A A«t 



-l+l/ag-x^ 



a 



— e 



r(l/a) A 



for A > and t > 0, on substituting x = z a to obtain the second equality. 
The final expression in (2.26) reveals a connection with the standard normal 
distribution corresponding to a = 2. Indeed, by the scaling property, it is no 
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restriction to assume that t 
a = 2 reads 



1 so that the final expression in (2.26) with 



(2.27) 



2 1 

^A 6 



x 2 



~* 2 dz : 



— — erfc(A) 
A 



for A > 0. By (1) in [9], page 265, one knows that 

1 



(2.28) 



L/ 



erfc(A)](x) 



7T 



and hence it follows that 



(2.29) 



L' 



-i 



r e A^ 

— — erfc(A) 
A 



x 2 /4 



-y 2 /* 



dy 



for x > 0. The density function / of Si obtained on the right-hand side of 
(2.28) and the distribution function F of Si given on the right-hand side 
of (2.29) coincide with the expressions obtained from the reflection principle 
Mi := max < t <i B t = law |£i|, which yields Si = law y/2M x = law >/2|Bi| = law 
I-B2I, where B = (Bt)t>o is a standard Brownian motion. 

3. By the scaling property, it is no restriction to assume that t = 1 in 
the sequel. Let F denote the distribution function of Si and let / denote 
the density function of Si. Note that the form of the Laplace transform on 
the right-hand side of (2.24) being combined with (2.26) implies that the 
density function exists (see the proof of Corollary 2 below for further detail) . 



Combining (2.19), (2.24) and (2.26), on using that 
since F(x) = Jq f(y) dy for x > 0, it follows that 



1 r 1 



Uf](X)](x)=F(x) 



(2.30) 



a 



r(i/a 

for x > 0. 

To simplify the notation, consider the equation 



dz 



(x) 



(2.31) 

for x > 0, where we set 
(2.32) 



g(x)=L^[G(A)](x) 



G(X) 



for A > 0. From (2.32), we see that G"(A) 

G'(X) 



(2.33) 



a 



- z " dz 

aA"- 1 G(A) 
1 



1 so that 



A 



A 2- 



G(A) + - = 



for A > 0. Since 1^ X [G' {X)\(x) = -xg(x), it follows that h^ 1 [C (X) / \]{x) = 
— Jq yg(y)dy for x > 0. Moreover, using (1) in [9], page 137, we see that 
L7 l [l/X 2 ~ a ](x) = l/(r(2 - a)^" 1 ) so that L^ 1 [G(A)/A 2 - a ](x) = (l/r(2 - 
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a)) Jo(g(y)/(x — y) a ~ 1 )dy for x > 0. Finally, we have L^" 1 [l/A](x) = 1 for 
x > 0. Hence, taking L^ 1 in (2.33), we find that 

(2.34) - f yg(y) dy - — f 9{y) dy + l = Q 

Jo T(2-a) Jo [x-y) a 1 

for x > 0. Noting that g(x) = (T(l/a)/a)f(x) for x > 0, we see that (2.34) 
reads 

a r f(y) , « 



f x ol f x fly) 

(2.35) J q yf(y) dy + f ^ y ^ j—^ dy ^ 



for x > and this is exactly equation (2.5). 

4. We will seek a solution to (2.35) of the form 



(2.36) f(x) = ]T a n 

n=0 



where (3 and 7 are constants to be determined. First, note that 
(2.37) f yf(y) dy = f^ a n f y^ +l dy = f^ 

JO ~~n JO ~r 



n=0 u n=0 



Pn + 7 + 2 



for x > 0. Second, by (3.191) in [10], page 333, and (6.2.2) in [1], page 258, 
we have 

(2.38) f y^ 1 (x - yf- 1 dy = x^ l B(^ v) = x^' 1 T ^ T ^ 



for (i > 0, v > and x > 0. It follows that 

f^yio (*- y )-^ = r(2^|>"/o (s-y)"- 1 ^ 

(2.39) 

_^ ar(/?n + 7 + l) ^ n+7 _ a+2 
- n ^ a "r(/?n + 7-a + 3) X 

for x > 0. Combining (2.35), (2.37) and (2.39), we find that (3 = a and 
7 = a — 2. Inserting (2.37) and (2.39) into (2.35) with these (3 and 7, we get 



(2.40) Y, ( anAn + a n+ iB n+1 )x a{n+1) + a Q B 

n=0 

where the constants A n and B n are defined by 



a 



r«A^ a 1 j o r(a(n + l)-l) 

(2.41) A n = — and B n = a- 



a(n + 1) r(cm + 1) 
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for n > 0. Prom (2.40) and (2.41), we find, by induction, that 
(2-42) a n = (-l) n " 1 1 ° a 

for n> 1, where a = l/(r(l/a)r(a - 1)). Inserting (2.42) into (2.36) with 
(3 = a and 7 = a — 2, and making use of (2.41), we obtain the series repre- 
sentation 

m 



1 



r(i/a) 



■l) n T(na + l)r((n - l)a + 1) • • • T(a + l)r(l) 
£ a 2 ^n[ T((n + Tja - l)T(na - 1) • • • T(2a - l)T(a - 1) 

x a(n+l)-2 



(2.43) x 53 



for x > 0. Using Stirling's formula r(ax + 6) ~ v / 2~7re _a:r (ax) aa;+fc_1 / 2 as 
x — > 00, where a > and 6 G R (cf. (6.1.39) in [1], page 257), it is readily- 
verified that \a n+ i/a n \ = 0(n 1_0 ) as n — > 00, whence the ratio test implies 
that the series in (2.43) converges absolutely for every x > and that / 
defined by (2.43) is a continuous function on (0, 00). Note, also, that only 
the leading term (l/(r(a — l)r(l/a))x Q-2 of the series is singular at zero, 
so we can integrate in (2.43) term by term over any finite interval in [0, 00). 
Finally, by induction over n > 0, using the fact that Y{x + 1) = xT{x) for 
x G R \ {0, — 1, — 2, . . .}, it is easily verified that the series representation 
(2.43) can be simplified to the form given in (2.8) above. 

5. We now show that / from (2.8) is a unique solution to the integral 
equation (2.5). For this, let us first note that since / satisfies (2.43) and hence 
solves (2.35), it follows that g = (T(l/a)/a)f solves (2.34). Assuming that g 
has a Laplace transform and taking the Laplace transform L on both sides 
of (2.34), we see that G = h[g] solves (2.33). The general solution to (2.33) 
is given by G(A) = ce A " + e A " f?° e~ za dz for A > 0, where c is a constant. 
In order to compute the Laplace transform of / defined in (2.43), we could 
attempt to interchange L and the sum and use the fact that L[x p ](A) = 
T(p + 1)/A P+1 for p > — 1 and A > 0. Using the ratio test, however, it is 
possible to verify that the resulting series diverges and therefore is not equal 
to L[/]. We note, however, that if we could show that / °° e~ Xx f(x) dx — ► 
as A — > 00, then we would have c = and (2.30) would imply that / from 
(2.43) is indeed the density function of Si, as claimed. 

Given this difficulty, we shall take a different tack and establish uniqueness 
of the solution to (2.35) in the class of functions that are locally integrable 
on [0,oo) and bounded on compact subsets of (0, 00) [these conditions are 
natural requirements so that the left-hand side of (2.35) makes sense]. Mul- 
tiplying both sides of (2.35) by (z — x) a ~ 2 and integrating the resulting 
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identity with respect to x from to z, we can use Fubini's theorem and 
(2.38) to obtain 

— [ Z y(z-yr- l f(y)dy + aT(a-l) f" f(y)dy 
(2.44) a ~ lJ ° 



a 



-z 



a-l 



{a- l)r(l/a)' 

for z > 0. Note that the interchange of the order of integration above is 
justified whenever / is locally integrable on [0, oo) and bounded on com- 
pact subsets of (0, oo). Differentiating this identity with respect to z and 
substituting x for z, we get 

1 f x y 11 

( 2 - 45 ) ~~r? TT / 7 ^/(j/)^ + /W = rn/ wv T\~2^a 

aT[a — 1) Jo [x — yY a T{l/a)T[a — 1) x 1 a 

for x > 0. This is a weakly singular Volterra integral equation of the second 
kind. Previous considerations show that both the function from (2.43) and 
the density function from (2.30) solve the equation (2.45). We note in passing 
that when / is the density function, we see from (2.45) that f{x) < [1/ (T(a — 
l)r (1 / a))](l / x 2 ~ a ) for all x > so that / is bounded on compact subsets 
of (0,oo). 

Denote by <f> the difference between the two solutions to (2.45). Then 

(2.46) a f ^2^</>(y) dy + 0(x) = 

Jo {x-y) z a 

for x > 0, where we set a = l/{aT{a — 1)). It follows from [12], Theorem 
7, page 35, that = if is locally square-integrable, but since the latter 
could not be the case (around zero), we give a direct proof of the former 
fact. For this, fix x\ > arbitrarily large and set £ = \<f>\. Letting 

(2-47) T£(x) = f t-^— dy, 



/o {x-yf 
we find, by induction using (2.46), that 
(2.48) f (z) < 6 n T n £(x) 

for x G (0, a;i] and n > 1, where b = ax\. An iterative calculation using Fu- 
bini's theorem and (2.38) shows that 

(x - y) 1 -^- 1 ) 

for x G (0, xi] and n > 1 with some constant c„, > 0. Choosing n > 1 large 
enough so that 1 — n(a — 1) < 0, combining (2.48) with (2.49) and applying 
a simple iteration procedure to the resulting inequality, we find that 



gnij,m-l 



( 2 - 5 °) £(x)< / £(y)dy 
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for x G (0, xi] and m > 1 with some constant c > 0. Since the right-hand side 
converges to zero as m — ► oo, it follows that £(a;) = for x € (0, si] and thus 
4>(x) = for all x > 0. This shows that the two solutions to (2.45) coincide 
on (0, oo). Hence, we can conclude that / from (2.8) is a unique solution 
to (2.5) in the class of functions which are locally integrable on [0, oo) and 
bounded on compact subsets of (0, oo). 

6. Note that the fractional differential equation (2.6) follows from (2.5) 
by differentiation so that (2.8) defines its solution satisfying the boundary 
condition (2.7). Now, suppose that / solves (2.6) and satisfies (2.7). Then 
(2.5) follows from (2.6) by integration [on using (2.38) with fi = a — 1], so / 
solves (2.35). Then proceeding as above, we find that / must be equal to the 
density function, as long as / is locally integrable on [0, oo) and bounded 
on compact subsets of (0, oo). This establishes the existence and uniqueness 
claim about (2.6) and (2.7) in the latter class of functions. The proof of the 
theorem is complete. □ 

Remark 1. The integral equation (2.5) is closely related to the (gener- 
alized) Abel equation of the first kind 

(2.51) £ (a + * ) f(y) dy = R(x) (0</3<l), 

which admits a closed-form solution expressed in terms of the Riemann- 
Liouville fractional derivative of R (of order 1 — j3). For more details, see 
[16] and the references therein. Note that the integral equation (2.5) is of 
the form 

(2.52) £ (ay + * ) f(y) dy = R{x) (0</3< 1), 

which may be viewed as being of the first order if the Abel equation (2.51) 
is viewed as being of the zeroth order. Note also that the equation (2.45) is 
the "second kind" analog of the equation (2.5). 

Remark 2. The results of Theorem 1 extend to the case when the Levy 
measure equals 

Q 

(2.53) v(dx) = -j^dx, 

where c > is a general constant. This can be derived using the scaling 
property of X. Letting, in this case, ft denote the density function of St = 
sup 0<s<t X s , we note for future reference that (2.8) extends as follows: 



(2.54) f t (x) = J2 



(cr(-a)t) n ~ 1 /«r(an - l)T(-n + 1 + 1/a 



-x 



an— 2 
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for x > and t > 0. Similarly, from (2.30), it is readily verified that 
(2.55) Ee~ XSt = e - Xx f t (x)dx = ——^e KtXa e~** dz 

for A > and t > 0, where we set n = cT(—a). Note, in particular, that (2.55) 
yields ES t = (a/T(l/a))(nt)^ a for t > 0. 



7. Further to the series representation given in (2.8) above, the next corol- 
lary presents an integral representation for the density function / of Si . Since 
this representation extends to t ^ 1 and c / 1/T(— a) by the scaling property 
of X, we will only focus on the case when t = 1 and c = 1/T(— a) in (2.53). 
We refer to [11], Theorem 1, page 422, and [6], Theorem 3, page 74, for more 
general integral representations in this context (with no obvious connection 
to the one given below). 



Corollary 2. Let X = {Xt)t>o be a stable Levy process of index a G 
(1,2) satisfying (2.1) and (2.2), and let S\ =sup < t <iA^ t denote its supre- 
mum over the time interval [0, 1]. Then the density function f of S\ is given 
by 



1 f c 

m = - / 



(2.56) 



t a cos(a7r/2) 



+ 



cos(t Q sin(a7r/2) + tx) 



t a ycos(ocir/2) 



T(l/a) Jo {t a -yY~ l / a 

x sin(y sin(a7r/2) + tx) dy 



dt 



for x > 0. 



Proof. Setting in the first equality and noting in the second equality 
that 

(2.57) H(X) = f°° e^ dy = - f* e""" dy 

J x a Jo 

for A > 0, it follows from (2.30) that 

(2-58) /(x) = _0_ L -i[ e A^ (A)](x) 

for x > 0. From the second equality in (2.57), one sees that H can be ana- 
lytically continued to the entire complex plane. The same fact is therefore 
true for A i— > e A ° H(\) so that the Laplace inversion formula is applicable in 
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(2.58) yielding 

r(l/a) 2m J-ioo V a 

2.59) 

= ^J_ r>>-(-!W2> _ a (\-mr dy ) e «* dt 



r(l/a) 27r J-oo V " 



o 



for x > 0. In the case t > 0, we have exp((zt) a ) = exp(t a (cos(a7r/2) + 
is'm(an/2))) and exp(— (ity) a ) = exp(— (ty) Q (cos(a7r/2) + i sin(a7r/2))). In 
the case t < 0, we have exp((zi) Q ) = exp((— i) Q (cos(a7r/2) — isin(a7r/2))) 
and exp(— (ity) a ) = exp(— (— ty) a (cos(air/2) — isin(air /2))). Inserting these 
expressions into (2.59), one can verify that the integral from — oo to +oo 
equals twice the integral from to +oo, which, in turn, can be reduced to 
the form given in (2.56) above. As this verification is somewhat lengthy, 
but still straightforward, further details will be omitted. This completes the 
proof. □ 

8. The next corollary describes asymptotic behavior of the law of Si at 
zero and infinity. Recall that f(x) ~ g(x) as x — > xq means that 
lima-^o f( x )h{x) = 1 for G [-oo,+oo]. 

Corollary 3. Let X = (Xt)t>o be a stable Levy process of index a £ 
(1,2) whose characteristic function is given by 

(2.60) Ee iXXt = exp(tJ™(e iXx - 1 - *Ax)^ 

for A G K and t > 0, where c > is a given and fixed constant. Let S\ = 
sup < t <i-Xt an d f denote the density function of Si. Then 

(2.61) fix) ~ ; — : ■ 1 ; ; , /' 2 US X j 0, 

V ; JK J (cr(-a)) 1 " 1 /«r(a- l)r(l/a) 

(2.62) f(x)^cx~ a - 1 as x]oo. 

Proof. The relation (2.61) follows directly from the explicit series rep- 
resentation (2.54). The relation (2.62) can be derived from the integral rep- 
resentation (2.56), as shown in [8]. □ 
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